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Using a sample of (225.3 ± 2.8) x 10 6 J/ip decays collected with the BESIII detector at BEPCII, 
searches for invisible decays of rj and rj' in J/ip — > <jyq and cprf are performed. Decays of <f> — > K + K~ 
are used to tag the r\ and rj decays. No signals above background are found for the invisible 
decays, and upper limits at the 90% confidence level are determined to be 2.58 x 10~ 4 for the ratio 

B(^invisiblc) g gg x 10 -2 f B(^->invi S iblc) 

B{r)— >77) o(j7'— 

PACS numbers: 13.25.Gv, 13.20.Jf, 14.40.Be 



I. INTRODUCTION 

Invisible decays of quarkonia to final state particles 
that are not observable in general purpose particle de- 
tectors offer a window into physics beyond the standard 
model (SM) [3-[3]- This is because, aside from neutrinos, 
the SM includes no other invisible particles that these 
states can decay into. Observation of the quarkonium 
states coupling to invisible final states can provide in- 
formation on new neutral particles such as photinos or 
neutralinos, very light gravitinos, or light dark matter 
matter candidates (LDM) [l|-|3|. These LDM particles 
may have the necessary density to account for the non- 
baryonic mass of the Universe [J] . 

It is in any case very interesting to search for such light 
invisible particles in collider experiments especially 
since LDM is a proposed interpretation for the origin of 
the 511 keV line from the galactic bulge observed by the 
SPI spectrometer on the International Gamma-Ray As- 
trophysics Lab (INTEGRAL) satellite @. Many searches 
for the invisible decays of 7r°, rj, ?/, J/ip and T(15) have 



been performed by previous experiments [7H10| . Invisible 
r\ (rj 1 ) decays may produce a pair of LDM particles in the 
MeV mass range, with expected rates depending on the 
spin and couplings of the particle mediating the decays. 
A spin-1 U boson [ll| would contribute only through its 
axial couplings to quarks. Invisible r\ (?/) decays could 
also occur from rj (rj 1 ) — > UU, with the U decaying in- 
visibly into neutrinos or LDM particles [2]. For spin-0 
boson exchanges tentative order of magnitude estimates 
for the branching fractions to a pair of LDM particles 
B(r){r}') -> XX) ~ 1-4 X 10" 4 (1.5 x 1CT 6 ) have been 
given in [12], using the WMAP [12] result on the relic 
density of the Universe. 

Using 58 x 10 6 J/ip events, the BESII experi- 
ment obtained a first upper limit of B(rj(rj') — > 
invisible)/^ (?/) 77) < 1.65 x 1CT 3 (6.69 x 
1CP 2 ), corresponding to an upper limit of B(r)(rf) — > 
invisible) < 6.5 x 1CT 4 (1.5 x 1CT 3 ) [8]. Complemen- 
tary to the BESII results, IceCube set B(rj — > v e ^ T v e ^) < 
6.1 x 10~ 4 [HI for V decays into SM neutrinos. 

In this paper, we present updated results of searches 
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for the invisible decays of 77 and 77'. The data sample used 
consists of (225.3 ± 2.8) x 10 6 J/tp events [H[ collected 
with the BESIII detector at the BEPCII collider [13]. 



II. THE BESIII EXPERIMENT AND MONTE 
CARLO SIMULATION 

BEPCII/BESIII is a major upgrade of the BESII 
experiment at the BEPC accelerator. The design peak 
luminosity of the double-ring e + e~ collider, BEPCII, is 
10 33 cm~ 2 s- 1 at a beam current of 0.93 A. The BESIII 
detector has a geometrical acceptance of 93% of 47r and 
consists of four main components: (1) a small-celled, 
helium-based main draft chamber (MDC) with 43 lay- 
ers, which provides measurements of ionization energy 
loss (dE/dx). The average single wire resolution is 135 
(Um, and the momentum resolution for charged particles 
with momenta of 1 GeV/c in a 1 T magnetic field is 0.5%; 
(2) an electromagnetic calorimeter (EMC) made of 6240 
Csl (Tl) crystals arranged in a cylindrical shape (barrel) 
plus two end-caps. For 1.0 GeV photons, the energy res- 
olution is 2.5% in the barrel and 5% in the end-caps, and 
the position resolution is 6 mm in the barrel and 9 mm 
in the end-caps; (3) a time-of-flight system (TOF) for 
particle identification (PID) composed of a barrel part 
made of two layers with 88 pieces of 5 cm thick, 2.4 m 
long plastic scintillators in each layer, and two end-caps 
with 96 fan-shaped, 5 cm thick, plastic scintillators in 
each end-cap. The time resolution is 80 ps in the barrel, 
and 110 ps in the end-caps, corresponding to a 2a K/ir 
separation for momenta up to about 1.0 GeV/c; (4) a 
muon chamber system made of 1000 m 2 of resistive-plate- 
chambers arranged in 9 layers in the barrel and 8 layers 
in the end-caps and incorporated in the return iron of 
the super-conducting magnet. The position resolution is 
about 2 cm. 

The optimization of the event selection and the estima- 
tion of physics backgrounds are performed using Monte 
Carlo (MC) simulated data s amp les. The GEANT4-based 
simulation software BOOST 18] includes the geometric 
and material description of the BESIII detectors, the de- 
tector response and digitization models, as well as the 
tracking of the detector running conditions and perfor- 
mance. The production of the Jjifi resonance is simu- 
lated by the MC event generator kkm c fl9j | ; the known 
decay modes are generated by EVTGEN [20| with branch- 
ing ratios set at PDG values [2l|, while the remaining 
unknown decay modes are modeled by lundcharm [22I ]. 

III. DATA ANALYSIS 
A. Analyses for 77 and rj' — ^ invisible 

In order to detect invisible 77 and 77' decays, we use 
J — > (fir] and (fir)'. These two-body decays provide a 
very simple event topology, in which the <fi candidates can 



be reconstructed easily and cleanly decaying into K + K~ . 
The reconstructed (fi particles can be used to tag 77 and 
7/ in order to allow a search for their invisible decays. In 
addition, both the <fi and 77(77') are given strong boosts in 
the J/tfi decay, so the directions of the 77 and 77' decays 
are well defined in the lab system and any decay products 
can be efficiently detected by the BESIII detector. The 
missing 77 and 77' can be searched for in the distribution 
of mass recoiling against the (fi candidate. 

Charged tracks in the BESIII detector are recon- 
structed using track- induced signals in the MDC. We 
select tracks that originate within ±10 cm of the in- 
teraction point (IP) in the beam direction and within 
1 cm in the plane perpendicular to the beam. The tracks 
must be within the MDC fiducial volume, | cos#| < 0.93 
(8 is the polar angle with respect to the e + beam di- 
rection). Candidate events are required to have only 
two charged tracks reconstructed with a net charge of 
zero. For each charged track, information from TOF 
and dE/dx are combined to calculate XpidM values. 
With the corresponding number of degree of freedom, 
we obtain probabilities, ProbpiD(i), for the hypothe- 
ses that a track is a pion, kaon, or proton, where i 
(i = ir/K/p) is the particle type. For both kaon candi- 
dates, we require ProbpiD(K) > ProbpiD^)- The mass 
recoiling against the (fi candidate, M^ oco11 , is calculated 
using the four-momentum of the incident beams in the 
lab frame {p^ ah = + ) , and constructing the 4- 
product (M^ CCDl1 ) 2 = (pi ab -pKKY(pub ~Vkk)p., where 
Pxk = P4, i s the sum of the four-momentum of the two 
charged kaons. The 77 and 77' signal regions in the MJ occul 
distribution are defined to be within 3cr of the known 
masses of 77 and 77' [2l|. Here, a is the detector resolution 
and is 17.8 (9.3) MeV/c 2 , which is determined from MC 
simulation, for J/tp —> ^77(77'). 

Electromagnetic showers are reconstructed from clus- 
ters of energy deposits in the EMC crystals. The shower 
energies are required to be greater than 25 MeV for the 
barrel region (| cos#| < 0.8) and 50 MeV for the end-cap 
region (0.86 < | cos6*| < 0.92). The showers in the tran- 
sition region between barrel and end-cap are required to 
have an energy greater than 100 MeV. Showers must be 
isolated from all charged tracks by more than 10°. 

We require that 77(77') — > invisible events have no 
charged tracks besides those of the <fi — > K + K~ candi- 
date. In addition, the number of EMC showers (A^hower)! 
that could be from a Kl or a photon, are required to 
be zero inside a cone of 1.0 rad around the recoil direc- 
tion against the <fi candidate. This requirement rejects 
most 77 and 77' decays into visible final states. It also 
eliminate most backgrounds from multibody decays of 
J jifi — > 0+anything. In order to ensure that 77 and 77' 
decay particles are inside the fiducial volume of the de- 
tector, the recoil direction against the (fi is required to 
be within the region | cos 6* rcco ii | < 0.7, where 6> roco ii is 
the polar angle of the recoil three-momentum of <fi can- 
didate. Figure [1] (a) shows the K + K~ invariant mass 
distribution after the above selection. A clear cfi peak 
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FIG. 1: (a) The vhkk distribution for candidate events in data. The arrows on the plot indicate the signal region of cj> 
candidates. Points with error bars are data; the (blue) histogram is expected background, (b) Recoil mass distribution against 
4> candidates, M™ cml , for events with 1.01 GeV/c 2 < tukk < 1-03 GeV/c 2 in (a). Points with error bars are data; the (blue) 
solid histogram is the sum of the expected backgrounds; the dashed histograms (with arbitrary scale) are signals of rj and rj 
invisible decays from MC simulations; the arrows on the plot indicate the signal regions of the rj and rj' — > invisible. 




FIG. 2: The M™ coil distribution with events around the rj' 
mass region. Points with error bars are data. The (black) 
solid curve shows the result of the fit to signal plus back- 
ground distributions, the (blue) dotted curve shows the back- 
ground shape from J/ip -> </>/o(980)(/o(980) -> K L K L ), 
the (blue) dashed curve shows the polynomial function for 
J/ip — > 4>KlKl background, and the (red) dotted-dash curve 
shows the signal yield. 



is seen. Figure [T] (b) shows the recoil mass against <p 
candidates for events with 1.01 GeV/c 2 < ttikk < 1-03 
GeV/c 2 , and there are no significant signals in the rj and 
rj' mass regions. 

We use MC simulated events to determine selection 
efficiencies for the signal channels and study possible 
backgrounds. The efficiencies are 36.0% and 36.1% for 
rj and rj' invisible decays, respectively. More than 20 
exclusive decay modes are generated with full MC sim- 
ulations in order to better understand the backgrounds. 
The sources of backgrounds are divided into two classes. 




M reooil (GeV/c 2) 

FIG. 3: The A^ ccoil distribution for the control sample J/ip -> 
4>rj' , rj' — > n + n~rj(r] — > 77) decay candidates. The solid curve 
shows the fit results. 



Class I: The background is from J/ip — > <pr)(r)'), where 
<p —> K + K~ and 77(77') decays into visible final states that 
are not detected by the EMC. The expected number of 
background events from this class is 0.18±0.02 (1.0±0.2) 
in the signal region for the 77(77') case. Class II: It is from 
J/ip decays to final states without 77(77') or without both 
77(77') and (p. For the 77 invisible decay, the dominant back- 
ground is from J/ij) — > jrj c , rj c — > K it^Kl, while for the 
77' case, the dominant background is from J/ip — > (pK^Kj, 
and J/tp -> </>/ (980), / (980) -> K L K L . The expected 
number of background events from class II is 0.8±0.2 and 
9.4 ± 1.7 in the signal regions for 77 and 77', respectively. 

After all selection criteria are applied, only one event 
(shown in Fig.[T](b)) survives in the 77 signal region where 
1.0 ± 0.2 background event is expected. An upper limit 
(UL) at the 90% confidence level (C.L.) of N^ L = 3.34 for 
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J/tp — > <f>r) ((f) — > K + K~ and 77 — > invisible) is obtained 
using the POLE ++ program |23l) with the Feldman- 
Cousins frequentist approach |24j . The information used 
to obtain the upper limit includes the number of observed 
events in the signal region, and the expected number of 
background events and their uncertainty. 

For the 77' case, an unbinned extended maximum like- 
lihood (ML) fit to the MJ ocoil distribution in the range 
0.8 GeV/c 2 < A/£ ocoil < 1.2 GeV/c 2 , as shown in Fig. |3 
is performed. The signal shape used in the fit, shown 
in Fig. [31 is obtained from a nearly background-free 
J/ip — > <p7]\ V — > 7r + 7r~?7, rj — > 77 sample. The 
purity of the sample is greater than 98.5%. The shape 
of the invisible signal peak in the Af^ eco11 distribution 
is fixed to the smoothed histograms of the J/ip — > <pr\' , 
rj — > 7r + 7r~?7, 77 — > 77 MC sample, and the signal yield 
is allowed to float. The shape of the dominant back- 
ground J/ip 0/o (980), / (980) K L K L is described 
by MC simulated data, in which the ,fn (980) line shape 
is parameterized with the Flatte form [25J 

/(m) = Ml+mi + z(gl P7r7r + g2p KK y (1) 

where Mf is the mass of the / (980), m is the effective 
mass, p is Lorentz invariant phase space (p = 2k /m, here, 
k refers to the n or K momentum in the rest frame of 
the resonance), and g\ and 32 are coupling-constants for 
the /o(980) resonance coupling to the irn and KK chan- 
nels, respectively. These parameters [Mf = 0.965±0.010 
GeV/c 2 , g\ = 0.165 ±0.018 (GeV/c 2 ) 2 and g\ = 0.695 ± 
0.075 (GeV/c 2 ) 2 ] have been determined in the analysis of 
J/tp -> 07T+7T- and <j)K+K- from BESII data daE^. 
In the ML fit, the dominant background shape (J/ip —> 
#0(980), / (980) -> K L K L ) is fixed to the MC simula- 
tions, and its yield (N* 5 ) is floated. The shape of the 
remaining background from J/ip — > (pK^Kj^ is modeled 
with a first order Chebychev polynomial whose slope and 
yield (N no ^_f o ) are floated in the fit to data. The signal 

yield, iVj- g = 2.3 ± 4.3, is consistent with zero observed 
events, and tne resulting ntted values ofiV^ kg andiVl k „ g f 

' ° to non-/o 

are 239 ± 28 and 37 ± 25, respectively, where the errors 
are statistical. We obtain an upper limit by integrating 
the normalized likelihood distribution over the positive 
values of the number of signal events. The upper limit 
at the 90% C.L. is N?! T = 10.1. 



The selection criteria for the charged tracks are the 
same as those for J/tp — > <jyq(r]'), 77(77') — > invisible. How- 
ever, at least two good photons are required. The events 
are kinematically fitted using energy and momentum 
conservation constraints (4C) under the J/ip — > KK^ 
hypothesis in order to obtain better mass resolution and 
suppress backgrounds further. We require the kinematic 
fit XK+K--yy to k e l ess than 90 (40) for the 77(77') case. If 
there are more than two photons, the fit is repeated using 
all permutations, and the combination with the best fit 
to KK^i^i is retained. 




$.35 ' 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 

r\ (GeV/c 2 ) 



FIG. 4: The (a) vtikk and (b) m 77 distributions with fit re- 
sults superimposed for J/ip — > 4>rj, <p —¥ K + K~, rj 77. 
Points with error bars are data. The (black) solid curves 
show the results of the fits to signal plus background, and the 
(black) dashed curves are for signal. In (a) , the (blue) dotted- 
dash curve shows non-0-peaking backgrounds, and the (red) 
short-dashed curve shows the non-77-peaking background. In 
(b), the (blue) dotted-dash curve shows non-r?-peaking back- 
grounds, and the (red) short-dashed curve shows the non-s- 
peaking background. 



B. Analyses for r\ and 77' —y 77 

The branching fraction of 77(77') — > 77 is also deter- 
mined in J/ip —¥ ^77(77'), in order to obtain the ratio of 
6(77(77') — > invisible) to 6(77(77') — > 77). The advantage of 
measuring )~^' nvisl ^ lc ) j s that the uncertainties due 

to the total number of J/ip events, tracking efficiency, 
PID, the number of the charged tracks, and the residual 
noise in the EMC cancel. 



The numbers of J/ip — > <jn](rf), 77(77') — > 77 events 
are obtained from an extended unbinned ML fit to the 
K + K~ versus 77 invariant mass distributions. The pro- 
jection of the fit on the rriKK (777-77) ax ^ s ^ s shown in 
Figs, da) andEJa) (Figs. Egb) and Mb)) for the 77 and 
77' cases, respectively. In the ML fits, we require that 
0.99 GeV/c 2 < m KK < 1.10 GeV/c 2 and 0.35 GeV/c 2 
< m 77 < 0.75 GeV/c 2 (0.75 GeV/c 2 < TO77 < 1.15 
GeV/c 2 ) for the 77(77') case. The signal shape for <f> is 
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FIG. 5: The (a) ttikk and (b) m 11 distributions with fit re- 
sults superimposed for J/rp — > cfyrf ', <f> — > K + K~, rj' — > 77. 
Points with error bars are data. The (black) solid curves 
show the results of the fits to signal plus background distribu- 
tions, and the (black) dashed curves are for signal. In (a), the 
(blue) dotted-dash curve shows non-c/i-peaking backgrounds, 
and the (red) short-dashed curve shows the non-r/-peaking 
background. In (b), the (blue) dotted-dash curve shows non- 
r/-peaking backgrounds, and the (red) short-dashed curve 
shows the non-^-peaking background. 



chev polynomial function (B(m 77 )). All parameters re- 
lated to the background shape are floated in the fit to 
data. The PDFs for signal and backgrounds are com- 
bined in the likelihood function £, denned as a function 
of the free parameters NJf 



^non-077 

Jv bkg ' 



AT non-4> 
iV bkg ' 



and N^: 



C = 



-(N 7 ! +N I, ° n " l ' T >+N I, ° D " , '+N' 1 ° a ~' 1 ) 



An 



N 



x \{[N^RBW{m KK ) x CB{m^ 

8=1 

+N^B{m KK ) x B(m^) 
+N™^B(m> KK ) x CB(m^) 
+N™pRBW(m< KK ) x B(m^)], 



(3) 



where 



-> 0?7, (p 



K+K- 

AT non-r) 
iV bk K 



V 



non-0 



and 



is the number of J/'ip 
-> 77 events, and A^™~ , ,» bkg 
are the numbers of the corresponding three kinds 
of backgrounds. The fixed parameter N is the total num- 
ber of selected events in the fit region, and m l KK (to 77 ) 
is the value of uikk ) for the ith event. We use the 
product of the PDFs, since we have verified that ttikk 



and 



are uncorrelated for each component. The nega- 



tive log-likelihood (— ln£) is then minimized with respect 
to the extracted yields. The resulting fitted signal and 
background yields are summarized in Table HI We also 
obtain the results for the rj' case by replacing 77 with 77' 
in Eq. ([3]). The fitted results for 77(77') — » 77 are shown in 
Fig- HI (Fig- E]). The detection efficiencies are determined 
with MC simulations to be 36.3% and 31.7% for 77 and 
77', respectively. 



modeled with a relativistic Breit-Wigner (RBW) func- 
tion [28| convoluted with a Gaussian function that rep- 
resents the detector resolution; the signal shape for 77(77') 
is described by a Crystal Ball {CB) function [2{|, and 
its parameters are floated. In the ML fits, the width 
of 4> is fixed at the PDG value, and its central mass 
value is floated. The backgrounds are divided into three 
categories: non-077(?7')-peaking background {i.e., J/rp — ¥ 
r fTr°K + K~ , in which one of the photons is missing); non- 
</>-peaking background {i.e., J/ip — > K + K~rj{r]')); and 
non-77(77')-peaking background {i.e., J/rp — > ^77 and 
</>7T 7r° ). The probability density functions (PDF) for 
non- </>-peaking background in the ttikk distribution is 
parameterized by [30] 

B{m KK ) = {m KK - 2m K ) a • r^—KK, ( 2 ) 

where a, b and c are free parameters, and ttlk is the nom- 
inal mass value of the charged kaon from the PDG [2l[ . 
The shape for the non-77(77')-peaking background in the 
m 77 distribution is modeled by a second-order Cheby- 



TABLE I: The fitted signal and background yields for J/ip — > 
4> r l{ r l l )j viv') 77> an d e 77( e 77 ) is its selection efficiency. 





Value 


Quantity 


V V 




13390 ± 136 400 ±25 


,>Tnon-^77 / A/ -non-c/>7/ \ 
JV bkg U v bkg ) 


2514 ± 64 1482 ± 46 


Ajnon-</>/ ii-non-^ 
Jv bkg l Jv bkg ) 


1132 ±70 10 ±15 


jl Tnon-77 / » rnon-77 \ 
iV bkg V iV bkg J 


313 ± 54 159 ± 26 


e v (e v ' ) 

c 77 V c 77/ 


36.3% 31.7% 



According to the results in Table HI the ratio of 
B{J/rp — > (j>r)) to B(J/rp — » <j>r)'), is found to be consistent 
with the known value [21| • The individual branching frac- 
tion is larger by 1.3[1.6)er with respect to the average 
value listed in Ref. for B{J/rp — > ,677(77')), while it is 
consistent with Ref. [81. 
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IV. SYSTEMATIC UNCERTAINTIES 

The contributions to the systematic error on the cal- 
culation of the ratios are summarized in Table [TTJ The 
uncertainty, due to the requirement of no neutral show- 
ers in the EMC inside the 1.0 rad cones around the recoil 
direction against the <p candidate, is estimated using the 
control sample of fully reconstructed J/tp — ¥ cpr), r\ — > 77 
events. The ratios of events with no extra photons to 
events without this requirement are obtained for both 
data and MC simulation. The difference 0.3% is taken 
as the systematic error for both the 77 and 77' cases. This 
study determines the difference of the noise in the EMC 
for MC simulation and data. The uncertainty due to the 
<j> mass window requirement is determined to be 1.5% by 
using the same control sample of J/tp — > cpr/, 77 — > 77 
events. 

For the 77' invisible decay, systematic errors in the ML 
fit originate from the limited number of events in the 
data sample and from uncertainties in the PDF param- 
eterizations. Since the signal shape is obtained from 
the J/tp — > <jyrf , T)' — > n + ir~r], 77—7-77 events in the 
data, the uncertainty due to the signal shape is negli- 
gible. The uncertainty due to the background shape is 
estimated by varying the PDF shape of the background 
in the ML fit. The shape of the dominant background 
J/tp -> </3/ (980), / o (980) ->■ K L K L is parameterized 
with the Flatte form in Eq. (fl}. To estimate the un- 
certainty, we change the central values of the param- 
eters used in the fit by one standard deviation of the 
measured values [26| . and find that the relative error on 
7/' — > invisible decay is 1.0%. The systematic uncertainty 
due to the choice of parameterization for the shape of 
the background from J/tp —> <PKlKl is estimated by 
varying the order of the polynomial in the fit; we find 
a relative change on the invisible signal yield of 2.9%, 
which is taken as the uncertainty due to the background 
model. 

The uncertainty in the determination of the number 
of observed J/tp —> <prj{rj'), <p — > K + K~, 77(77') — s- 77 
events is also estimated. The systematic error due to 
photon detection is determined to be 1% for each photon 
[3l| . The uncertainty due to the 4C fit is estimated to be 
0.4%(0.8%) for the 77(77') case using the control sample 
J/tp — > ir°K + K~. In the fit to the <f> mass distribu- 
tion, the mass resolution is fixed to the MC simulation; 
the level of possible discrepancy is determined with a 
smearing Gaussian, for which a non-zero a would repre- 
sent a MC-data difference in the mass resolution. The 
uncertainty associated with a difference determined in 
this way is 0.1% (1.0%) for the 77(77') case. The system- 
atic uncertainty due to the choice of parameterization for 
the shape of the non-<^77(?7')-peaking background is esti- 
mated by varying the order of the polynomial in the fit; 
we find the relative changes on the 77(77') signal yield of 
0.1% (0.6%), which is taken as the uncertainty due to 
the background shapes. The total systematic errors cr^ ys 
and o\?, s on the ratio are 2.5% and 4.1% for 77 and 77', as 



summarized in Table [TT1 

TABLE II: Summary of errors. The first four lines are rel- 
ative systematic errors for J/tp —¥ ^77(77'), 77(77') — > invisible. 
The next four lines are relative systematic errors for J/tp — > 
'Pvi 7 )'): viv') ~* 77- The second line from the bottom is the 
relative statistical error of NZj y (NIt„). 







Sys. 


error (%) 


Source of uncertainties 




V 


v' 


Requirement on JV a hower 




0.3 


0.3 


<p mass window 




1.5 


1.5 


Background shape of J/tp — > 


0/o (980) 




1.0 


Background shape of J/tp — > 


<j>K L K L 




2.9 


4C fit for 77(77 ) — > 77 




0.4 


0.8 


Photon detection 




2.0 


2.0 


Signal shapes for 77(77') —5- 77 




0.1 


1.0 


Background shape for 77(77') - 


-» 77 


0.1 


0.6 


Total systematic errors 




2.5 


4.1 


Statistical error of N^(N!^) 




1.0 


6.0 


Total errors 




2.7 


7.4 



V. RESULTS 

The upper limit at the 90% confidence level on the 
ratio of 6(77 — > invisible) to 6(77 — > 77) is calculated with 

6(77 -> invisible) ^ Nujjey 1 ^ 
B{r) -> 77) Nljj/e^j 1 - a n ' 

where N^j L is the 90% upper limit of the number of ob- 
served events for J/ip — > <fyq, <p — > K + K~ , 77 — > invisible 
decay, e v is the MC determined efficiency for the signal 
channel, NJj is the number of events for the J/tp — > cprj, 
(p — > K + K~, i] — > 77, elj 7 is the MC determined 
efficiency, and is the total error for the 77 case from Ta- 
ble |TTJ The upper limit on the ratio of B{rf — > invisible) 
to 73(77' — > 77) is obtained similarly. 

Thus, the upper limit of 2.58 x 10" 4 (2.39 x 10" 2 ) on 
the ratio of ,6(77(77') — > invisible) and 6(77(77') — » 77) is 
obtained at the 90% confidence level. 



VI. CONCLUSION 

In summary, the invisible decays of 77 and rf are 
searched for in the two-body decays J/tp — > (prj and <pr]' 
using (225.3 ± 2.8) x 10 6 J/tp decays collected with the 
BESIII detector. We find no signal above background 
for the invisible decays of 77 and 77' and obtain upper 
limits at the 90% C.L. of 2.58 x 10~ 4 and 2.39 x 10~ 2 
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for 



BOq— ^-invisible) 



and 



B(W— ^invisible) 
B(?7'->77) 



, respectively. Using 



the branching fraction values of 77 and 77' — » 77 from the 
PDG [2 1[| , we determine the invisible decay rates to be 
B(r] —> invisible) < 1.0 x 10~ 4 and B{j]' — > invisible) < 
5.3 x 1CT 4 at the 90% confidence level. 

Our limits are improved by factors of 6 and 3 compared 
to the previous ones obtained at BESII 8(, the 77' limit 
being almost 2 times better than the recent one from the 
CLEO-c experiment [Hj]. The limit for 77 — > invisible 
is smaller than a possible estimate [l2| for the 77 — > XX 
decay to a pair of MeV dark matter particles, no such 
decays, however, being expected for a spin-1 U medi- 
ator with vector couplings to quarks (or axial ones in 
the limit of small m x ). These limits constrain the de- 
cays 77 (77') — > UU with each U decaying invisibly (into 
neutrinos or LDM), with branching fraction B mv . The 
resulting 77 (77') limits on the U couplings to quarks are 
improved by ~ 1.6 and 1.3 as compared to those ob- 
tained in [2] from the BESII limits @, and now read 
V/2 + ft < 3 ^- 2 /VB~ and |/ s | < 4 lO- 2 /VB~, 
respectively, where f u , fd and f s are effective couplings 
of U boson to light quarks 0j. 
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